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Supersymmetric neutralinos are one of the most promising candidates for the dark matter in the 
Universe. If they exist, they shonld make np some fraction of the Milky Way halo. We investigate 
the fluxes of positrons expected at the Earth from neutralino annihilation in the halo. Positron 
propagation is treated in a diffusion model including energy loss. The positron source function 
includes contributions from both continuum and monochromatic positrons. We End that, for a 
“canonical” halo model and propagation parameters, the fluxes are generally too low to be visible. 
Given the large uncertainties in both propagation and halo structure, it is however possible to obtain 
observable fluxes. We also investigate the shapes of the positron spectra, including fits to a feature 
indicated by the results of the Heat experiment. 
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I. INTRODUCTION 

It is well known that a large fraction of the mass of the Universe has only been observed by its gravitational effects. 
The nature of this dark matter is unknown. Standard Big Bang Nucleosynthesis bounds the baryon density of the 
Universe, but these results fall far short of the known amount of dark matter. Thus more exotic forms of matter are 
sought, and weakly interacting massive particles (WIMPs) are one of the most promising cold dark matter candidates. 

Perhaps the best WIMP candidate is the neutralino (for a review, see Ref. [I]), which arises in supersymmetric 
models as a linear combination of the superpartners of the neutral gauge and Higgs bosons. We choose the neutralino 
as our WIMP candidate, though we allow its properties to vary over a generous sample of supersymmetric models. 

It may be possible to detect neutralinos (or WIMPs) in the galactic halo by the products of their mutual annihi¬ 
lations, e.g. by searching for gamma rays [2], antiprotons [3] and positrons [4-6] coming from the Milky Way halo. 
In this paper we consider neutralino annihilation in the halo into positrons (both continuum and monochromatic). 
Compared to earlier studies [4-6], we will use a true diffusion model instead of a leaky-box model. This has the 
advantage of allowing us to use more realistic models for the diffusion constant while still permitting an analytic 
solution. The numerical implementation is very fast, making detailed scans over the huge MSSM parameter space 
feasible. We also include all two-body final states of neutralino annihilation (at tree level) using detailed Monte Carlo 
simulations for the decay and/or hadronization of the annihilation products. We will compare our results mainly with 
those of Kamionkowski and Turner [6] (hereafter KT). 

We will also compare our predictions of the positron fluxes to observations. The fluxes are typically small compared 
to the expected background, but given the large uncertainties involved, we identify models that can have a significant 
effect on the observed positron spectrum. We will show that a large class of models improves the fit of observations 
to the predicted background. Finally, we investigate the positron spectra at higher energies than have been observed, 
and identify features that may be detectable by future experiments such as Ams [7]. 


II. SET OF SUPERSYMMETRIC MODELS 

We work in the Minimal Supersymmetric Standard Model (MSSM). In general, the MSSM has many free parame¬ 
ters, but with some reasonable assumptions we can reduce the number of parameters to the Higgsino mass parameter 
fi, the gaugino mass parameter M 2 , the ratio of the Higgs vacuum expectation values tan/3, the mass of the CP-odd 
Higgs boson 771 , 4 , the scalar mass parameter rrio and the trilinear soft SUSY-breaking parameters Ab and At for the 
third generation. In particular, we don’t impose any restrictions from supergravity other than gaugino mass unifi¬ 
cation. We have made some scans without the GUT relation for the gaugino mass parameters Mi and M 2 . This 
mainly has the effect of allowing lower neutralino masses to escape the LEP bounds and is not very important for 
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TABLE I. The ranges of parameter values used in our scans of the MSSM parameter space. Note that several special scans 
aimed at interesting regions of the parameter space has been performed. In total we have generated approximately 103000 
models that are not excluded by accelerator searches. 


this study. Hence, the GUT relation for Mi and M 2 is kept throughout this paper. For a more detailed definition of 
the parameters and a full set of Feynman rules, see Refs. [8,9]. 

The lightest stable supersymmetric particle is in most models the lightest neutralino, which is a superposition of 
the superpartners of the gauge and Higgs fields, 

X° = NiiB + Ni2W^ + Ni 3H° + NuHl ( 1 ) 

It is convenient to define the gaugino fraction of the lightest neutralino, 

Zg = \Nii\^ + \Ni2\^. (2) 

For the masses of the neutralinos and charginos we use the one-loop corrections as given in Ref. [10] and for the Higgs 
boson masses we use the leading log two-loop radiative corrections, calculated within the effective potential approach 
given in Ref. [11]. 

We make extensive scans of the model parameter space, some general and some specialized to interesting regions. 
In total we make 18 different scans of the parameter space. The scans are done randomly and are mostly distributed 
logarithmically in the mass parameters and in tan /3. For some scans the logarithmic scan in fj, and M 2 has been 
replaced by a logarithmic scan in the more physical parameters and Zg/{1 — Zg) where is the neutralino mass. 
Combining all the scans, the overall limits of the seven MSSM parameters we use are given in Table I. 

We check each model to see if it is excluded by the most recent accelerator constraints, of which the most important 
ones are the LEP bounds [12] on the lightest chargino mass, 

f 91 GeV , I'm + — 'm^o \ > 4 GeV 

m ^ ' Xi vi I ('o'l 

^1 \ 85 GeV , otherwise 

and on the lightest Higgs boson mass nifjo (which range from 72.2-88.0 GeV depending on sin(/3 — a) with a being 
the Higgs mixing angle) and the constrains from & ^ 57 [13]. For each allowed model we calculate the relic density 
of neutralinos where is the density in units of the critical density and h is the present Hubble constant 

in units of 100 km s“^ Mpc“^. We use the formalism of Ref. [14] for resonant annihilations, threshold effects, and 
finite widths of unstable particles and we include all two-body tree-level annihilation channels of neutralinos. We also 
include the so-called coannihilation processes in the relic density calculation according to the analysis of Edsjo and 
Gondolo [ 8 ]. 

Present observations favor h = 0.6±0.1, and a total matter density Um = 0.3±0.1, of which baryons may contribute 
0.02 to 0.08 [15]. Not to be overly restrictive, we accept in the range from 0.025 to 1 as cosmologically interesting. 

The lower bound is somewhat arbitrary as there may be several different components of non-baryonic dark matter, but 

we demand that neutralinos are at least as abundant as required to make up the dark halos of galaxies. In principle, 
neutralinos with < 0.025 would still be relic particles, but only making up a small fraction of the dark matter 
of the Universe. We will only consider models with < 0.025 when discussing the dependence of the signal on 

III. POSITRON SPECTRA FROM NEUTRALINO ANNIHILATION 

When neutralinos annihilate in the galactic halo they produce quarks, leptons, gauge bosons, Higgs bosons and 
gluons. When these particles decay and/or hadronize, they will give rise to positrons either directly or from decaying 
mesons in hadron jets. We thus expect to get both monochromatic positrons (at an energy of 'm^) from direct 
annihilation into e+e“ and continuum positrons from the other annihilation channels. In general the branching ratio 
for annihilation directly into e“''e“ is rather small, but for some classes of models one can obtain a large enough 
branching ratio for the line to be observable. 
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For continuum positrons, we simulate the decay and/or hadronization with the Lund Monte Carlo Pythia 6.115 
[16]. We perform the simulation for a set of 18 different neutralino masses from 10 to 5000 GeV and for the “basic” 
annihilation channels consisting of the heavy quarks, leptons and gauge bosons. The annihilation channels containing 
Higgs bosons are then easily taken into account by allowing the Higgs bosons to decay in flight and averaging the 
produced positron flux over the decay angles. For any given MSSM model, the positron spectrum is then given by 


dfj) 

de 


d(j) 


de 


cont. 



^ Bp 

F^e+e“ 


d(j) 

de 


+ Bg+g- S{e 

F 


rnx), 


( 4 ) 


in units of e“''/annihilation where e = Eg+ /(I GeV), rhx = rrix/iX GeV), Bp is the branching ratio into a given final 
state F and d(l)lde\p is the spectrum of positrons from annihilation channel F. We include all two-body final states 
(except the three lightest quarks which are completely negligible) at tree level and the Zy [17] and gg [18] final states 
which arise at one-loop level. 


IV. PROPAGATION MODEL 


A. Propagation and the interstellar flux 


We consider a standard diffusion model for the propagation of positrons in the galaxy. Charged particles move 
under the influence of the galactic magnetic field. With the energies we are concerned with, the magnetic gyroradii 
of the particles are quite small. However, the magnetic field is tangled, and even with small gyroradii, particles can 
jump to nearby field lines which will drastically alter their courses. This entire process can be modeled as a random 
walk, which can be described by a diffusion equation. 

Positron propagation is complicated by the fact that light particles lose energy quickly due to inverse Compton and 
synchrotron processes. Diffuse starlight and the Cosmic Microwave Background (CMB) both contribute appreciably 
to the energy loss rate of high energy electrons and positrons via inverse Compton scattering. Electrons and positrons 
also lose energy by synchrotron radiation as they spiral around the galactic magnetic field lines. 

Our detailed treatment of positron diffusion follows. First we define a dimensionless energy variable e = FjiX GeV). 
The standard diffusion-loss equation for the space density of cosmic rays per unit energy, dn/de, is given by 


d dn 
dt de 


K{e,x)V 


dn 

de 


de 


b{e, x) 


dn 

de 


Q{e,x), 


( 5 ) 


where K is the diffusion constant, b is the energy loss rate and Q is the source term. We will consider only steady 
state solutions, setting the left hand side of Eq. (5) to zero. 

We assume that the diffusion constant K is constant in space throughout a “diffusion zone”, but it may vary with 
energy. At energies above a few GeV, we can represent the diffusion constant as a power law in energy [19], 

K{e) = Ao£“ « 3 X 102^£° ®cm2 s'^ (6) 


However, at energies below about 3 GeV, there is a cutoff in the diffusion constant that can be modeled as 

K{e) = ATo [C -b £“] « 3 X 10^^ [3° ® -b e° ®] cm^ s'^ (7) 

We will in the following consider both of these models for the diffusion constant, but we will focus on the second 
expression. We will refer to Eq. (6) as model A and Eq. (7) as model B. The function 6(e) represents the (time) rate 
of energy loss. We allow energy loss via synchrotron emission and inverse Compton scattering. The rms magnetic field 
in the diffusion zone is about 3 /xG, an energy density of about 0.2 eV cm“^. We allow inverse Compton scattering 
on both the cosmic microwave background and diffuse starlight, which have energy densities of 0.3 and 0.6 eV cm“^ 
respectively. These two processes combined give an energy loss rate [20] 

6(£)e± = —e^ « 10“^®e^ s“\ (8) 

Tp 

where we have neglected the space dependence of the energy loss rate. The function Q is the source of positrons in 
units of cm“^ s“^. 

The diffusion zone is a slab of thickness 2L. We will fix L to be 3 kpc, which fits observations of the cosmic ray flux 
[19]. We impose free escape boundary conditions, namely that the cosmic ray density drops to zero on the surfaces of 
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the slab, which we let be the planes z = ±L. We neglect the radial boundary usually considered in diffusion models. 
This is justified when the sources of cosmic rays are nearer than the boundary, as is usually the case with galactic 
sources. We will see that the positron flux at Earth, especially at higher energies, mostly originates within a few kpc 
and hence this approximation is well justified in our case. 

We will now consider the diffusion-loss equation, Eq. (5), in the steady-state. We first make a change of variables, 
u = 1/e, and we write the diffusion constant as a function of u, with K{E) = KQh{u), yielding 

1 d dvi _ dvt [- 9,/ \i — 1^/ / \ —#\ 

AM a;* = (9) 

where te = 10^® s is the energy-loss time introduced in Eq. (8). We now suppose that a < 1, and we make a further 
change of variables v = f h{u)du. We find the following transformations for models A and B, 

VA = - v}~°' and vb = Cu + - - (10) 

I — a I — a 

These transformations need to be inverted. The first one can be done analytically, 

u = [{1 - a)vA]~ , ( 11 ) 

though the second one requires a numerical inversion. We notice that the function of u in the source term is just the 
Jacobian of the full transformation, 


1 

u^h{u) 

Again, in the first case this can be evaluated. 



2-q 

w{va) = [(1 - a)vA] , 


( 12 ) 


(13) 


and the second requires the implicit solution of u{v). We put F{v) = dn/du(v) and arrive at 

^F(v, x) = KoteV‘^F{v, x) - TEwiv)Q{e{v),x). (14) 

Eq. (14) is equivalent to the inhomogeneous heat equation. The space variables are exactly analogous, while the 
variable v takes the place of time in the heat equation. We also notice that the source function is multiplied by a 
weight factor —w{v). 

The interpretation of the “time” variable v is worth a short discussion. The variable v increases monotonically with 
decreasing energy. Since this model does not include any reacceleration of positrons, the variable v will only increase. 
This is precisely the situation of the time variable in the standard heat equation. 

In addition to the boundary conditions previously stated, we must also impose an “initial” condition. It is clear 
that the initial condition should he F(v = 0,x) = 0. This implies that the positron density is zero at ?; = 0, which is 
infinite energy. 

With these things in mind, we can write the solution as an integral over the source multiplied by a Green’s function. 
We transform back to the original energy variable, remembering that Eq. (14) is an equation for dn/du and not dn/de. 
The solution can be written 

e) r 

dv'w(v') / d^x' G2l[v{£^) — v',x — x']Q(£(v'),x') (15) 

where the function G 2 L is the Green’s function for Eq. (14) on the slab with thickness 2L. By our construction, 
w(v)dv = —de and making this variable change frees us from the cumbersome Jacobian factor and from implicitly 
inverting the transformation v(e). The only place v{e) remains is in the argument to the Green’s function. The formal 
solution is thus written 


dn 

de 


= TeS 


-2 


dn 

de 


TEe ^ J J G 2 L[v{e) — v{e'),x — x']Q{e',x'). 


The free-space Green’s function for Eq. (14) is given by 


(16) 
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Ghee{v-v',x-x') = [iTTKoTE{v-v')\ exp 6»(v - v')■ (17) 

V iK^TEiv - v') J 

We require a Green’s function that vanishes on the planes z = ±L, the boundaries of our diffusion region. We can 
use a set of image charges, 


y'n = y\ 


4 = 2Ln+(-l)"z', 


(18) 


to find the Green’s function, 


OO 

G2l{v-v',X-x')= ^ {-l)^Gfree{v-v',X-x'^). 

n= —OO 


(19) 


It is simple to verify that this Green’s function vanishes on the boundaries. 

Now we choose an appropriate source function Q. To simplify the calculation, we will assume that the source 
is uniform in z and depends only on the cylindrical coordinate r. Let the number density of WIMPs in the halo 
be nog{rsph), where rgph is the spherical coordinate and no is the number density in the center of the galaxy. The 
annihilation rate, and thus the source strength, is proportional to the square of the density. We define 


f{r) 


1 

TL 


l-L 


dz g (^sph)- 


We will mainly use the isothermal halo profile. 


5iso(r,z) = 


'sph 


( 20 ) 


( 21 ) 


with a =5 kpc and a local halo density of 0.3 GeV/cm^. We will refer to this is our “canonical” halo model. We will 
also compare with the Navarro, Frenk and White profile [21], 


5nfw(?’, 2 ;) =---3 (22) 

Fsph (rsph + a) 

with a = 25 kpc. Finally, the source function is 

Q{e,r) =nlf{r){av)tot^, (23) 

where (crv)tot is the annihilation rate and d(j)/de is the positron spectrum in a single annihilation. This contains 
annihilations both to monochromatic and continuum positrons. Note that if < 0.025, we need to rescale the 
source function by (n^/i^/0.025)^. We will only include models with < 0.025 (with the rescaled source function) 
in Fig. 3 (b). 

We now insert this source into the Green’s function integral, Eq. (16), 


dn 

de 


LLIl 2 / \ 

— = nQ{av)totTE£ 


de' [inKoTEivie) - u(e'))] X7 X (“1) 


- 3/2 d(j) 


00 


de' 

n— — oo 



^27r 


dr'r'f{r') / dO' exp 



2rr' cos6>' + (z-(-l)"z' 
4iFoT£(u(£) - v{e')) 



(24) 


The z' integration will simply yield error functions and the 9' integration yields 27r/o, where /q is the modified Bessel 
function of the first kind. We are left with the r' integral, which we put in the function X, 


I[r, z, Au] 


1 


AKqTe^V 


E Ee^'f 


n= —cxD ± 


(-l)”L + 2Ln±z 
^/AKqteAv 


dr' r'f{r')io 


/ 2rr' \ 


4KoTeAv ) 


(25) 
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FIG. 1. (a) The effective diffusion time tb as a function of An. When doubling the diffusion constant K 27 , positrons remain 
trapped for a shorter time before escaping the galaxy. At low energy losses (small An) this of no importance since it is the 
energy loss time that is important in that region. When doubling the energy loss time rie we increase the flux at low energy 
losses since the positrons don’t lose energy as quickly. At higher energy losses, the flux goes down since the positrons have 
now had more time to escape. Increasing L to 5 kpc, we see that diffusion becomes effective at a larger An, as expected. The 
positrons require more time to escape the larger diffusion region. The 10% difference at low An is because the function /(r) 
depends on L. In (b) we show the effective diffusion time td{xE, E) versus the injection energy E for a:=0.50, 0.25 and 0.10. In 
(c) we show the effective diffusion time td{xE,E) as a function of x for different injection energies. In (b) and (c) propagation 
model B is used. 


with Io{x) = Io{x)e Note that when x ^ 1, Io{x) « I\/2ttx. When AKqTe^v <C this sum and integral is 
easy to perform, since we have a simple gaussian and all of the error functions cancel. The result is 

X [r, z, A?; < L^j{^K qTe)] = /(r). (26) 

For clarity, we define an energy dependent diffusion time, 

td[Au] = TD[e, e' > e]= teX [?;(£) - n(£')] . (27) 

This encodes all of the information about the spatial extent of the source function and the energy dependence of the 
diffusion constant. What this function represents is the lifetime against energy loss and diffusion of positrons at a 
given energy which were emitted at a higher energy. There are two regimes to this timescale. When An is small, the 
energy lost is small, thus the particles have not been propagating for long, and have not had any chance to escape. 
At large An, much energy has been lost, thus the particles have been propagating for a long time and have had a 
chance to escape the diffusion zone. In the former regime, it is the energy loss that dominates the propagation, and 
in the latter, it is the diffusion. Diffusion occurs in both regimes, but in the former, it occurs over a small enough 
region that the boundary can be neglected. At intermediate Au, there is a slight rise in the function td, before it 
drops precipitously. This is because the effective amount of source material increases with increasing Au. As more 
energy is allowed to be lost, the positrons can come from farther away. This allows us to sample the higher density 
inner regions of the galaxy. Only when the boundary can be probed does te begin to drop. 

In Fig. 1 (a) we show te as a function of An for the position of the solar system, r = 8.5 kpc and z = 0 kpc, 
for various combinations of K27 = Kq/IO'^'^ cm^ s“^, tiq = te/W^^ s, and L. At low Av, te is independent of K27 
and linear in te, as this is the regime where energy loss dominates. At high An, te decreases with increasing K 27 , 
as expected for the diffusion regime. Increasing te implies that particles have been propagating longer, and have 
had more chance to escape. Increasing L to 5 kpc increases the signal at high Av since it is more difficult to escape 
the diffusion zone. At low Av, the 10% discrepancy is due to the fact that f{r) depends on L. The small size of 
the discrepancy indicates that using a source function uniform in z is adequate. In Fig. 1 (b) we show the function 
te{xE,E) versus the injection energy E for j:=0.50, 0.25 and 0.10. In Fig. 1 (c) we show te{xE,E) versus x for 
some different injection energies E. In both Fig. 1 (b) and (c), propagation model B is used. 

The positron spectrum is now given by 

^ = nl{av)tot£~^ J de' ^^e [e, e '], (28) 

giving the total positron spectrum 
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(29) 


dn 

de 


{crv)tot£ BuneTD[£,mx] + 




td[£,£'] >0 {rhx - e) ■ 


Remembering that this is an expression for the number density of positrons, the flux is given by 

d^ Pc dn c dn 
de Att de Att de ’ 


(30) 


where /3c is the velocity of a positron of energy e. For the energies we are interested in, /3c ~ c is a very well justified 
approximation. 

We find that it is numerically advantageous to perform the integral in the variable v. Replacing the Jacobian factor, 
we find 


dn 

de 


= nl{av)tot£ BuneTD[£,mx] + 


/ dv w{v j 


de 


[e{v')]TD[e,e{v')] U (m^ - e), 


(31) 


The function td = can be tabulated once. The integrand is smooth, and is thus easily computed for a range of 
values, equally spaced in In Ac. Likewise the sum on error functions converges rapidly, and for the range of Ac values 
we are concerned with, need not be taken past n = ±10. 


B. Solar modulation 


The formalism developed above gives the interstellar flux of positrons from neutralino annihilations. There is a 
further complication in that interactions with the solar wind and magnetosphere alter the spectrum. These effects 
are referred to as solar modulation. This can be neglected at high energies, but at energies below about 10 GeV, the 
effects of solar modulation become important. We will not consider solar modulation, since its effects can be removed 
by considering the positron fraction, e+/(e+ ± e“), instead of the absolute positron fluxes. 

However, another complication is the possibility that the solar modulation is charge-sign dependent as indicated 
in Ref. [22]. Following their treatment we can write the positron fraction at Earth in a solar cycle and A~ 
respectively as 


P{R + l)-f 




/e — 


i?(2/-l) 

p{R+l)-fR 

2/-1 


(32) 

(33) 


where / is the interstellar positron fraction and R is the ratio of the total flux of positrons and electrons in a solar 
cycle to that in an A~ cycle. The value R can be measured, and is approximately given by [22] 


R{e) = max [min (0.45 ±0.17 lne,l),0.18]. 


(34) 


We notice that in the absence of charge-sign dependence, the positron fraction is unaffected by solar modulation 
effects. As we will see later, the charge-sign dependent solar modulation worsens the agreement between the mea¬ 
surements and the expected background and given the large uncertainties involved, we will not use these expressions 
when predicting the positron fraction from neutralino annihilation. 


C. Positron and electron backgrounds 

When we want to compare our predictions with experiments, we must consider the expected background of positrons. 
Since most experiments measure the positron fraction and not the absolute positron flux, we also need to consider 
the electron background. We will use the most recent calculation by Moskalenko and Strong [25] (hereafter MS). For 
their model 08-005 without reacceleration, we have parameterized their calculated primary electron and secondary 
electron and positron fluxes as follows: 
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FIG. 2. The positron fraction for the background is shown both without and with charge-sign dependent solar modulation 
[22]. The normalization has been kept free and is fitted to the Heat 94-1-95 data [23]. Also shown is the Caprice 94 data [24]. 
The backgrounds are parameterizations to the results by MS [25]. 
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GeV ^ cm ^ s ^ sr ^ 


(35) 

(36) 

(37) 


where we as before have used e = E((l GeV). These parameterizations agree with their curves to within 10-15% for 
the whole intervals given in Ref. [25] (approximately 0.001-1000 GeV for primary electrons and 0.01-100 GeV for the 
secondary electrons and positrons). The parameterizations also have the same asymptotic slopes as the calculated 
fluxes at both the low and high energy end and extrapolations should be acceptable if not too far out of the regions 
given above. These predictions also agree roughly with the absolute flux measurements by the Heat experiment [26]. 
The experimental error bars are however smaller on the measurements of the positron fraction and in Fig. 2 we have 
made a fit to the Heat 94-b95 measurements [23] of the positron fraction from the background alone keeping only the 
normalization of the positron flux free. We also show the best fits with charge-dependent solar modulation included. 
The normalization factors for the best fits are kb = 1.11 and kt = 0.609 without charge-dependent solar modulation 
and with A'^ modulation respectively (bl+ is the correct state for the Heat measurement). The corresponding reduced 
x2 are 3.14 and 10.8, i.e. not very good fits. We also show the Caprice 94 measurements of the positron fraction 
[24]. 

We clearly see that the background is too steep at low energies and including the charge-sign dependent solar 
modulation only worsens the fit. For models with reacceleration, the fit is worse still [25]. We also note that there 
is an indication in the Heat data of an excess at 6-50 GeV. However, one should keep in mind that (a) the 
uncertainties below 10 GeV are large due to our poor knowledge of the solar modulation effects, (b) it is possible 
to make a better fit by hardening the interstellar nucleon spectrum [25] (even though this might be in conflict with 
antiproton measurements above 3 GeV [27]) and (c) it is possible to make a better fit by hardening the primary 
electron spectrum (even though this might be in conflict with direct electron measurements at higher energies [25]). 
We will not consider these issues here. Instead we will assume the background as given by the parameterizations above 
without charge-dependent solar modulation. We will then compare the positron background with our predictions and 
we also make a simultaneous fit of the positron fraction background signal from neutralino annihilations having only 
the normalization of the background and signal positron fluxes as free parameters. 
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FIG. 3. The flux of positrons from neutralino annihilation in the halo versus (a) the neutralino mass and (b) the neutralino 
relic density Q,^h?. In (a) only models with 0.025 < < 1 have been included and in (b) models with < 0.025 have 

been rescaled by {Q,^h^/0.025)^. The Heat 94 data in this energy interval is shown together with the background prediction 
by MS [25] (as given by our parameterization, Eq. (37)). 


To conclude the discussion of the backgrounds, there is clearly room (or even an indication) of an excess of positrons 
at intermediate energies that might be due to neutralino annihilations. However, the uncertainties are large and other 
explanations are feasible. We will entertain the possibility that the excess is due to neutralino annihilation in the 
next section. 


V. RESULTS 

We are now ready the present our results and compare with the expected background and experimental measure¬ 
ments. We will first consider the absolute fluxes and investigate the dependence on the propagation parameters. 
We will use our propagation model B with the isothermal sphere as our default model. We will also show some 
typical spectra. In the next subsection we will then discuss the predicted positron fractions and compare with the 
experiments, in particular the excess at 6-50 GeV indicated by Heat data [23]. 


Example 

number 

Model 

number 



Parameters and units 



X properties 


7^ 

GeV 

Mi 

GeV 

tan/3 

1 

ruA 

GeV 

mo 

GeV 

Ahimo 

1 

Atimo 

1 

mv 

GeV 

z, 

1 

1 

1 

JE28..005683 

-852.3 

-670.1 

13.1 

664.0 

1940.6 

-1.56 

-1.60 

335.7 

0.9951 

0.029 

2 

JEsp3.000426 

2319.9 

4969.9 

40.6 

575.1 

2806.6 

2.0 

0.27 

2313.3 

0.028 

0.16 

3 

JE29__034479 

-1644.3 

-202.2 

54.1 

181.7 

2830.9 

2.57 

2.50 

101.6 

0.99927 

0.042 

4 

JE23__000159 

-221.8 

-324.5 

1.01 

792.2 

2998.7 

1.71 

0.67 

130.3 

0.660 

0.058 


TABLE II. Example models. The parameter values are given together with the neutralino masses, gaugino fractions and 
relic densities. 
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Optimal Positron Energy, E (GeV) 


FIG. 4. The flux of positrons, <E>g+i?® versus the optimal positron energy, E, for which the ratio of the signal flux to the 
background flux is highest. The positron background from MS [25] (as given by our parameterization, Eq. (37)) is shown as 
the dashed line. Above 100 GeV, the background is an extrapolation of the results by MS. 


A. Absolute fluxes and spectra 

As a representative energy we will choose the average flux in the energy range 8.9-14.8 GeV which corresponds to 
one Heat bin (with average energy of about 11 GeV). In this energy range, we should be fairly unaffected by solar 
modulation effects. In Fig. 3 we show the absolute positron flux for propagation model B versus the neutralino mass 
and versus the relic density We also compare with the Heat 1994 measurement of the absolute flux [26]. We 

clearly see that with our canonical halo model we find no models that give fluxes as high as those measured. However, 
one should keep in mind that we probably have an overall uncertainty of about an order of magnitude coming from 
uncertainties regarding the propagation. We also have uncertainties coming from our lack of knowledge of the halo 
structure. For example, the halo could be clumpy which could easily increase the signal by orders of magnitude [28,29]. 
The local halo density is also uncertain by at least a factor of 2 (which corresponds to a factor of 4 uncertainty in 
the positron flux). In Fig. 3 (a) we see that once we are above the W mass, the spread of the predictions is much 
smaller. We also see that for heavy neutralinos, Higgsinos and mixed neutralinos typically give higher fluxes than 
the gaugino-like neutralinos. This is because the annihilation cross section to gauge bosons is very suppressed for 
gaugino-like neutralinos. For Higgsinos, the cross section to gauge bosons is also suppressed, but usually dominates 
anyway. 

In Fig. 3 (b) we see that the highest fluxes are approximately proportional to This is because is 

approximately proportional to {uv)~^ where (cm) is the thermally averaged annihilation cross section. The highest 
fluxes for a given are obtained when the annihilation proceeds mainly to channels giving high positron fluxes 

and the velocity dependence of the cross section is small (i.e. the thermally averaged cross section is close to the 
annihilation cross section at rest). In that case, it is essentially the same cross section that determines both and 
the positron fluxes and hence the strong (anti)correlation. When < 0.025 we rescale the flux by (r2^/i^/0.025)^ 

and hence the fluxes for low are proportional to Note that in all other figures, we will only show the 

models where 0.025 < < 1. 

In Fig. 4 we show the positron flux versus the optimal positron energy, where we define as the energy 

at which ((i4>/(ii?)signai/(d4>/(ii?)background is maximal. We also show the positron background prediction by MS [25]. 
We see that in many cases it is advantageous to look at higher energies than present experiments do. Again, we see 
that for heavier neutralinos, the mixed and Higgsino-like neutralinos give the highest fluxes. 

In Fig. 5 we show an example of a spectrum when the neutralino mainly annihilates into W'^W~. This model is 
example 4 in Table H. We choose this model as an example to explore the dependence of parameters and propagation 
models since it has two nice bumps from both primary as well as secondary and higher decays/hadronizations of the 
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FIG. 5. An example of a light model annihilating mainly into In (a), the differences between our models A and B 

and the KT model with energy dependent escape time [6] are shown. In (b) the effects (for model B) of increasing K, te, a 
and L are given and in (c), the differences between different halo profiles are shown. Note that the scale is different in (c). 


W bosons. In (a) we compare our model A and B with the propagation model of KT [6] with energy dependent 
escape time. Compared to our model B, the KT propagation model assumes (i) that the diffusion constant coefficient 
a = 1, (ii) that there is no low energy cutoff of the diffusion constant (i.e. as our model A), (hi) that the energy-loss 
time te is higher, (iv) that the diffusion constant Kq is higher and (v) that the disk is homogeneous. They also use 
a leaky-box model instead of a diffusion model. 

In Fig. 5 (a), we see that our model A gives higher fluxes than model B below about 10 GeV. This is easy to 
understand since model A has a lower diffusion constant at low energies and hence the positrons stay around longer 
and we get higher fluxes. On the other hand, the KT propagation gives a sharper peak from the primary decay of W~^ 
to e+. This is, as we will see in figure (b) and (c), due to the different values of Kg, te and a. At low energies, the 
KT model has about the same shape as our model A, since neither te nor the exact value of a are of big importance 
there. The different normalizations at low energies are due to different values of Kg. 

In Fig. 5 (b) we show what happens in model B if we change the value of Kq, te, a. or L. Increasing Kg reduces the 
flux at lower energies (where diffusion is important), whereas increasing te increases the flux at higher energies (where 
energy-loss is more important than escape through diffusion). We actually saw this behavior already in Fig. 1 (a) 
where td was shown versus Av for different values of Kg and r. Increasing a essentially tilts the spectrum clockwise 
at lower energies where diffusion is important. Hence, the dip at intermediate energies is more pronounced. Increasing 
L has more or less the same effect as lowering Kg, except at higher energies, where we see a decrease in the flux. This 
is merely an artifact of the averaging over z in Eq. (20). When L increases, the average at a given r decreases. We 
saw this behavior earlier in Fig. 1 (a). 

In Fig. 5 (c) we show the dependence on the actual halo profile. We expect to mainly see the positrons that are 
created within a few kpc of us, but as we go lower in energy (and hence higher in energy loss), the visible region 
expands. In (c) we show the difference between the isothermal sphere and a homogeneous disk. The isothermal sphere 
gives higher fluxes at low energies, i.e. we see the galactic core at lower energies. If we use a steeper halo profile like 
the Navarro, Frenk and White profile [21], the flux at lower energies goes up even more. The dependence on the 
precise form of the halo profile is never very large though. It is mostly the propagation parameters that determine 
the shape of the spectrum. 

To conclude, the differences between our models and the KT propagation model are mainly a matter of different 
parameters. 

In Fig. 6 we show some other typical spectra (examples 1-3 in Table II). In (a) we show a spectrum when gauge 
bosons don’t dominate in the final state. In this case, we typically get a big wide bump without any special features. 
In (b), we show a spectrum for a heavy neutralino annihilating into gauge bosons. As in Fig. 5, we see two bumps 
from the decay of the Ws. The right one comes from direct decay to and occurs at approximately The left 

one comes from secondary decays of Ws, r“*'(/r+) ^ e’*' and positrons from quark jets, W~^ ^ ^ 7r+ ^ e+. 

This bump occurs at approximately m^/60. Exactly where the bumps occur depend on the propagation model used. 
In (c), we show a spectrum of a light neutralino where we have artificially boosted the branching ratio for direct 
annihilation into e“'"e“ to 1%. Typically this branching ratio is below 10“^, but if we e.g. have a large mass splitting 
of the selectrons, i.e. mgj ^ and ^ we can find mostly gaugino-like, but also some mixed models with 
a branching ratio Bg+g- > 0.01. In this case, the line from monochromatic positrons might be visible and not too 
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FIG. 6. Example of positron spectra versus the positron energy. In (a) we show a typical spectrum for medium-heavy 
neutralinos, in (b) we show a spectrum from a heavy neutralino when gauge bosons dominate the annihilation final states and 
in (c) a spectrum when the branching ratio for annihilation into e'''e~ has been increased to 0.01. In all spectra we show them 
as calculated with our two models and with the energy-dependent model by KT [6] 


widened by propagation. We find that it is only when the branching ratio into e“'"e“ directly is reasonably high 
0.001) that the line is visible. We remark that if selectrons and smuons mix, the selectron degeneracy must be 
less than a few percent to avoid constraints from /x —> ey [30,1], though our model is allowed if the mixing between 
selectrons and smuons is negligible. 


B. Positron fractions 

Instead of comparing with the absolute fluxes, we will now compare with the measurement of the positron fraction by 
the Heat experiment [23]. As we saw in Fig. 2, the standard background predictions fail to reproduce the observations 
for all energies. Note, however, that the fit can be made better by adopting either a harder electron primary spectrum 
or a harder interstellar nucelon spectrum [25,27]. A harder electron spectrum seems to be in disagreement with 
high-energy electron measurements [25] and a harder nucleon spectrum might be in contradiction with antiproton 
measurements above 3 GeV [27]. In light of this, it is therefore interesting to investigate if the possible excess at 6-50 
GeV in the Heat data can be explained by neutralino annihilations in the halo. Note, however, the the uncertainties 
in the solar modulation (even though the charge-dependent solar modulation seems to worsen the fit) are large and 
not well understood. 

We will in the following entertain the possibility that the standard prediction of the background is correct (to 
within a normalization) and make a simultaneous fit of the normalization of the background and the signal. We 
will not include the charge-dependent solar modulation. A study of the Heat data in terms of primary sources 
from WIMP annihilations has previously been performed in Ref. [31] where the KT results (both with and without 
energy-dependent escape time) were used. 

We will compare with the most recent Heat measurement of the positron fraction (1994 and 1995 combined data) 
[23]. The error bars are smaller and the data cleaner for the positron fraction measurements. The Heat Collaboration 
gives their results in 9 bins from 1.0 to 50.0 GeV. 

As we saw in Fig. 3 (a), the positron fluxes are typically an order of magniture or more smaller than the Heat 
measurements, and we find that we need to boost the signal by a boost factor, kg, between 6 and 10^*^ to obtain a 
good fit. ks = 10^° is hardly realistic, but kg up to 100-1000 might be acceptable given the uncertainties in both 
propagation and halo structure (the halo could e.g. be clumpy [28,29]). 

In Fig. 7 we show some examples of our positron fractions compared with the Heat measurements. The examples 
are the same as those shown earlier for the absolute fluxes, i.e. examples 1-4 in Table H. The signals have been 
boosted by the boost factors given by the best fit. Almost all of our models have the general feature of increasing 
the flux at intermediate energies as shown in (a). For some models, where the annihilation predominantly occurs 
to gauge bosons, we also expect a peak at roughly m^/2 as shown in (b). The most pronounced feature, however, 
is the positron line. Even though smeared by propagation and energy loss, we still expect a sharp peak at if 
the branching ratio for annihilating into monochromatic positrons is high enough. However, the branching ratio is 
typically less than 10“^ and the feature is buried in the background. If the branching ratio is higher, the feature 
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FIG. 7. Example of positron fractions versus the positron energy. We show both the signal and background for the best 
simultaneous fit and also the best fit of the background alone. In (a) we show a typical spectrum for medium-heavy neutralinos, 
in (b) we show a spectrum for a heavy neutralino when gauge bosons dominate the annihilation final states, in (c) a spectrum 
when the branching ratio for annihilation into e''"e“ has been increased to 0.01 and in (d) a spectrum from a medium-heavy 
neutralino when Ws dominate. 


would be very clear as seen in (c). This might happen if there is a large mass splitting of the selectrons as discussed 
earlier. In (d) we show a light model annihilating mainly into gauge bosons and as seen, the bump from primary 
decay IT+ ^ e+ can be made to fit the excess quite nicely. 

In Ref. [31] models were found where the low-energy bump from gauge bosons can fit the excess at 6-50 GeV (as in 
our Fig. 7 (c)) without having to go to > 1 TeV neutralinos. They get better fits for lower masses than we do because 
they use the KT propagation model [6] which doesn’t shift features down in energy as much as our propagation does. 
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FIG. 8. The flux of positrons versus that of (a) antiprotons and (b) continuum gammas. The experimental limits from Bess 
95 [32] and high-altitude Egret [33] measurements are shown. 


C. Comparison with other signals 

We have seen that the positron fluxes are typically at least factors of a few too low to be observable. However, with 
only a small amount of boosting of the signal (coming either from different propagation models and/or clumping of 
the dark matter) we obtain positron fluxes that can either roughly fit the indication of an excess seen in the Heat 
data or produce features possible to observe with future experiments. We must be concerned with other signals that 
might be boosted at a similar level and that might be in conflict with current observations. 

In Fig. 8 we show the absolute flux of positrons versus the flux of (a) antiprotons [34,29] and (b) continuum gammas 
[29]. The antiproton measurement of the Bess collaboration [32] is shown as well as the limit of the high galactic 
altitude diffuse gamma emission as measured by Egret [33]. We see that especially the antiprotons are expected 
to give about the same or better constraints on neutralino dark matter than the positrons. For the models with 
reasonable boost factors (^100-1000) the antiproton flux is usually about a factor of 1-10 closer to the experimental 
bound than the positrons. We might therefore worry about getting antiproton fluxes that are too large when boosting 
the positron fluxes. However, even though a factor of 10 might seem to be high, we have large uncertainties involved 
in the propagation and solar modulation. For the positrons we also have large uncertainties from the energy loss that 
don’t enter in the antiproton propagation. Changes in the diffusion constant also affect the positron and antiproton 
fluxes differently. Increasing Kq, the positron fluxes at higher energies are mainly unchanged (see Fig. 5 (b)), whereas 
the antiproton fluxes decrease. We also tend to sample a larger volume for antiprotons than we do for positrons and 
different halo profiles and clumpiniess hence enter differently for antiprotons and positrons. A factor of 10 difference 
is thus probably within the uncertainties. 

For our examples and boost factors in Fig. 7, the corresponding maximal boost factors that would not violate the 
Bess upper limit [32] of 2.3 x 10“® GeV“^ cm“^ s“^ sr“^ on the antiproton flux would be kp = 9.8, 3.5 x 10^, 3.2 
and 10.7. Compare this with the boost factors from our best fits to the Heat data, kg = 11. 7, 1.1 X 10^, 19.0 and 
54.6. The last two examples. Fig. 7 (c) and (d) have boost factors which contradicts the antiproton measurements by 
a factor of 5-6 but as explained above, the uncertainties are large and we prefer to keep these kind of models anyway. 

The continuum gammas are very well correlated with the positron fluxes and of the same order of magnitude below 
present limits. The continuum gamma fluxes don’t suffer from large propagation uncertainties, but do depend more 
on the actual halo profile. 

To conclude, it is intriguing that both the positron, antiproton and continuum gamma fluxes for our canonical 
halo and propagation models end up reasonably close to the observed fluxes. The correlation is also quite good, even 
though different uncertainties enter in all three different predictions. 

Other signals, like direct detection of neutralinos and indirect detection from high-energy neutrinos from annihilation 
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in the Sun/Earth are not affected by dumpiness to the same extent as the ones from halo annihilation. Both direct 
searches and neutrino telescopes have started cutting into the MSSM parameter space considered here. However, the 
limits are not completely watertight and we choose to keep all models even though some of our models might be 
excluded with typical assumptions about the halo density and velocity dispersion. 


VI. CONCLUSIONS 

We conclude that for our canonical halo and propagation model, the positron fluxes from neutralino annihilation in 
the halo are lower than those experimentally measured. However, the predictions can easily be orders of magnitude 
too low due to our lack of knowledge of the structure of the halo (the halo could e.g. be clumpy). If we allow for this 
uncertainty we can easily obtain fluxes of the same order of magnitude as those measured. The shape of the signal 
spectrum is also different than the background and in most cases allows for a better fit of the excess at 6-50 GeV 
indicated by the Heat measurements than the background alone does. 

If positrons from neutralinos make up a substantial fraction of the measured positron fluxes below 50 GeV, we 
usually have some features to search for at higher energies, especially if the annihilation occurs mainly to gauge 
bosons where we expect a clear bump at m^l2. When the annihilation doesn’t go to gauge bosons, we only excpect a 
slight break in the spectrum at approximately mp^/2. This can be hard to find, unless we have detailed measurements 
around the break. If the branching ratio for annihilation into monochromatic positrons is higher than > 0.001, as can 
happen if there is a pronounced mass splitting between selectrons, there is a very pronounced feature to search for at 
a positron energy of my^. 

We also note that our propagation models don’t preserve features as well as that of KT [6]. This is primarily due 
to the fact that KT uses an escape time that implies a larger diffusion constant. Features are also shifted to lower 
energies with our propagation models. 

To conclude, there is a possibility to obtain measureable fluxes of positrons from neutralino annihilation in the 
Milky Way halo and they can easily be made to fit the excess indicated by Heat data. Future experiments will 
determine if there are features in the positron spectrum from neutralino annihilation. 
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